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Abstract 

Recent approaches to quark confinement are reviewed, with an emphasis on their connection 
to renormalization group methods. Basic concepts related to confinement are introduced: the 
string tension, Wilson loops and Polyakov lines, string breaking, string tension scaling laws, center 
symmetry breaking, and the deconfinement transition at non-zero temperature. Current topics 
discussed include confinement on R 3 x S 1 , the real-space renormalization group, the functional 
renormalization group, and the Schwinger-Dyson equation approach to confinement. 
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I. INTRODUCTION 



Quarks are the fermionic constituents of the hadrons, the strongly interacting particles 
such as the proton and the neutron. The problem of quark confinement, at its simplest, is to 
explain why quarks are not observed as physical states in high-energy scattering, but remain 
confined inside hadrons. The problem has been apparent for decades, andpre-dates quantum 
chromodynamics (QCD), the modern theory of the strong interactions l|-|3j]. Lattice gauge 
theory |4|, |5| has given us a physical basis for quark confinement: the force between quarks 
is independent of distance at large separation, and an infinite amount of energy would be 
required to isolate a single quark. The force constant at infinite separation is the string 
tension a. While lattice simulations have determined the relation of a to other observables 
of QCD with high precision, we are still far from a satisfactory understanding of quark 
confinement. We have neither a firm grasp of the mechanisms of a quark confinement nor 
the ability to calculate analytically even simple aspects of the physics of confinement. 

Recent work in the study of confinement suggests that progress is now being made, and 
there are several complementary approaches under development. All of these approaches 
have common features: they use the renormalization group in an essential way, and are 
strongly conditioned by lattice gauge theory results. This review is based on talks given 
at the INT Workshop "New applications of the renormalization group in nuclear, particle, 
and condensed matter physics", held February 22-26 2010, and the emphasis is on current 
research rather than a comprehensive treatment of each approach. The review is organized as 
follows: Section II introduces the basic concepts involved in the study of quark confinement: 
the string tension, its measurement, string breaking, the finite-temperature deconfinement 
transition, center symmetry, and string tension scaling laws. Section III discusses recent 
progress in understanding confinement on R 3 x S 1 Section IV describes recent work on a 
rigorous proof of confinement using real-space renormalization group methods. Section V 
discusses two related approaches to confinement based on the infrared behavior of quark 
and gluon propagators, the first using the functional renormalization group and the second 
using the Schwinger-Dyson equations. A final section attempts to summarize our current 
understanding and prospects for future progress. 
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II. FUNDAMENTALS OF CONFINEMENT 



The fundamental theory of the strong interactions is quantum chromodynamics (QCD), 
a quantum field theory describing the interactions of quarks and gluons. QCD is a gauge 
theory: Gluons are massless spin-1 bosons whose interactions are tightly constrained by the 
requirement of local gauge invariance in a manner similar to the case of QED (quantum 
electrodynamics) and photons. In QED, the symmetry group is U(l), the group of local 
phase transformations on charged fields: <f)(x) — > e ia ^(f)(x). In QCD, the corresponding 
symmetry group is 577(3), and the symmetry acts on the three color charges carried by 
quarks, usually labeled as red, blue and green. There is a profound difference between 
QED and QCD: the gauge group U(l) of QED is Abelian, but the gauge group SU(3) of 
QCD is non-Abelian. Physically, this implies that the photon has no electric charge but 
gluons do carry color charge. The four- vector field of the photon is written as A^x) where 
\i = 1,2,3,4, but an SU(N) gauge field is written as A a {x), where the index a runs over 
the N 2 — 1 one charges of the non-Abelian gauge field. While the photon is a free field with 
a trivial renormalization group beta function, SU (3) gauge theory, even without quarks, is 
a theory with highly non-trivial interactions. The beta function of an SU (N) gauge theory 
without quarks is 

to lowest order in perturbation theory. Here g is the strong-interaction coupling constant. 
It is actually a running coupling constant <?(//), depending on the momentum-space renor- 
malization scale fi, evolving according to 

^ = 0{9M)- (2) 

The negative sign of /3 (g) indicates that SU (N) gauge theories are asymptotically free: 
There is an ultraviolet fixed point at g = 0, and the interactions become weak at high 
energies. There is a corresponding infrared fixed point at g — oo, indicating that the running 
coupling g(fi) becomes large at low energies. This implies that the low-energy properties 
of QCD, unlike QED, cannot be analyzed via perturbation theory. The most effective tool 
we have for exploring the non-perturbative aspects of QCD is lattice gauge theory, in which 
continuous space-time is replaced by discrete lattice, with a fundamental cut-off scale set 
by the lattice spacing a. The renormalization group plays a crucial role in the taking of 
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the continuum limit, where a — > 0. In gauge theories, the gauge field has four-components 
but two physical polarization states. Two of the four components of the field correspond to 
redundant degrees of freedom, and must be eliminated by a choice of gauge, such as Coulomb 
gauge or Landau gauge, for continuum field theory methods to be applied. Lattice gauge 
theory, on the other hand, maintains exact gauge invariance: no gauge fixing is required, 
although it may be applied to lattice gauge theories. Physical predictions of theory, such 
as hadron masses, must of course be independent of the choice of gauge. A great deal of 
the technical problems in analytic treatments of non-perturbative aspects of gauge theories 
are ultimately associated with achieving gauge invariance and correct renormalization group 
scaling behavior simultaneously. 

The static potential between two heavy quarks can be determined in a gauge-invariant 
way from the expectation value of the Wilson loop. This is a non-local operator associated 
with a closed curve C in space-time parametrized as (r) where r can be taken to be the 
unit interval and x^ (1) = x^ (0). The Wilson loop is defined as 



where V indicates path-ordering of the gauge fields (x) along the path C. Taking W [C] to 
be an element of the abstract gauge group G, the basic observables are Tr R W [C], where the 
trace is taken over an irreducible representation R of G. TtrW [C] has a physical interpre- 
tation of the non-Abelian phase factor associated with a heavy particle in the representation 
R moving adiabatically around the closed loop C. Typically, we are interested in rectangular 
Wilson loops of with sides L and T. The loop can be associated with a process in which a 
particle-antiparticle pair are created at one time, move to a separation L, propagate forward 
in time for an interval T, and then annihilate. For T ^> L, we have 



where V R (L) is the heavy quark- ant iquark potential for particles in the representation R. 
For large distances, the potential can grow no faster than linearly with L, and that linear 
growth defines the string tension <j R for the representation R: 



as L — > oo. The O (1) correction term represents a so-called perimeter law contribution, 
which is not physical; there are also l/L corrections which are of physical interest. The 
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(3) 



{Tr R W[C}) ~exp[-ML) T] 



(4) 



V (L) ^ a R L + O (1) 



(5) 



statement that quarks are confined is the statement that the string tension op in the fun- 
damental representation of SU (3) gauge theory is non-zero. A characteristic feature of pure 
gauge theories in four dimensions is dimensional transmutation: the classical action is scale- 
invariant, but the introduction of a mass scale \x is necessary in order to define the running 
coupling constant g{p). It follows directly from the renormalization group that any physical 
mass in four-dimensional pure gauge theories must be proportional to a renormalization 
group invariant mass, usually written simply as A. This implies that or = crA 2 , where cr 
is a pure number. 

Paradoxically, string tensions are generally determined in so-called pure gauge theories, 
where only the gauge fields are dynamical; quarks and particles exist only as static classical 
charges. This is done because of string breaking , also known as charge screening. Consider 
a theory with dynamical particles of mass m in the fundamental representation F of the 
gauge group. Then energy obtained by separating a pair of static sources in the fundmental 
representation is opL\ when that energy becomes on the order of 2m, it will become ener- 
getically favorable to produce a pair of dynamical particles at an energy cost 2m, and no 
string tension will be seen for L > 2m/ a p. 

If one or more directions in space-time are compact, the string tension may also be 
determined using the Polyakov loop P, also known as the Wilson line. The Polyakov loop 
is essentially a Wilson loop that uses a compact direction in space-time to close the curve 
using a topologically non-trivial path in space time. The typical use for the Polyakov loop is 
for gauge theories at finite temperature, where space-time is R 3 x S 1 . The partition function 
being given by Z = Tr e~^ H , the circumference of S* 1 is given by the inverse temperature 
(3 = 1/T. In this case, we write 

f/3 



P (x) — V exp 



f 

i / dx^Ai (x) 
Jo 



(6) 



and string tensions may be determined from a two-point function 

{Tr R P (x) Tr R P^ (y)) ~ e -"*l*H7l (7) 

a behavior that assumes that the one-point function [TtrP (x)) = 0. The Polyakov loop 
one-point function {TtrP (x)) can be interpreted as the Boltzman factor exp (—(3Fr), where 
Fr is the free energy required to add a static particle in the representation R to the system. 
Of course, {TtrP (x)) = implies that F R — oo , which is thus a fundamental criterion 
determining whether particles in the representation R are confined. 



One of the most important concepts in our understanding of confinement is the role of 
center symmetry. The center of a Lie group is the set of all elements that commute with 
every other element. For SU(N), this is Z(N). Although the Z(N) symmetry of SU(N) 
gauge theories can be understood from the continuum theory, it is easier to understand from 
a lattice point of view. A lattice gauge theory associates link variable (x) with each lattice 
site x and direction /x. The link variable is considered to be the path-ordered exponential 
of the gauge field from x to x + fx: (x) = exp [iA^ (x)] . Consider a center symmetry 
transformation on all the links in a given direction on a fixed hyperplane perpendicular to 
the direction. The standard example from SU(N) gauge theories at finite temperature is 
U4 (x,t) — > zU^ (x,t) for all x and fixed t, with z G Z(N). Because lattice actions such as 
the Wilson action consist of sums of small Wilson loops, they are invariant under this global 
symmetry. However, the Polyakov loop transforms as P (x) — > zP (x), and more generally 



where k R is an integer in the set {0, 1, N — 1} and is known as the iV-ality of the rep- 
resentation R. If k R 7^ 0, then unbroken global Z(N) symmetry implies (Tr R P(x)) = 0. 
Thus global Z(N) symmetry defines the confining phase of a gauge theory. For pure gauge 
theories at non-zero temperature, the deconfinment phase transition is associated with the 
loss of Z(N) symmetry at the critical point T^. Below that point (TrpP (x)) = but above 



Notice that the case of zero iV-ality representations is special within this framework: 
there is no requirement from Z(N) symmetry that these representations are confined. This 
includes the adjoint representation, the representation of the gauge particles. However, 
lattice simulation indicate that (Tr R P (x)) is very small for these representations in the 
confined phase. Although screening by gauge particles must dominate at large distances, 
these zero iV-ality representations have well-defined string tensions at intermediate distances 
scales, e.g., on the order of a few fermi for SU(3), behaving in a manner very similar to 
representations with non-zero iV-ality {(], Q] • 

It is known from lattice simulations that each representation apparently has its own 
distinct string tension at intermediate distance scales. To a good approximation, the scaling 
behavior observed in lattice simulations is consistent with Casimir scaling 



Tr R P (x) ->■ z- 



kR Tr R P (x) 



(8) 



T d , (Tr R P(x))^0. 




(9) 
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where Cr is the quadratic Casimir operator for the representation R. This behavior is 
exact in two-dimensional pure gauge theories. It is generally believed that at sufficiently 
large distance scales, screening by gauge particles will lead to a single string tension for all 
representations within a given iV-ality class, i.e., with the same value of k. Labeling the 
lightest value of <7r within an iV-ality class by a^, Casimir scaling predicts 

a k k(N- k) 



ai N-l 

However, the observed behavior is also consistent with so-called sine-Law scaling 

Qk _ sin (f) 
ox sin h) 

and lattice results have not yet been definitive [8]. 



(10) 



(11) 



III. CONFINEMENT ON i? 3 x S 1 

In the last few years, it has proven possible to construct four- dimensional gauge the- 
ories for which confinement may be reliably demonstrated using semiclassical methods 
(9), [lo|]. These models combine Z(N) symmetry, the effective potential for P, instantons, 
and monopoles into a satisfying picture of confinement for a special class of models. All 
of the models in this class have one or more small compact directions. Models with an 
R 3 x S 1 topology have been most investigated, and discussion here will be restricted to this 
class. For recent work in other compactifications, see fill ]. The use of one or more compact 
directions will cause the running coupling constant of an asymptotically free gauge theory 
to be small, so that semiclassical methods are reliable. For example, if the circumference L 
of S 1 is small, i.e., L <C A -1 , then g(L) <C 1. However, this leads to an immediate problem: 
generally speaking, one or more small compact directions lead to breaking of Z(N) symme- 
try in those directions. This has been clear for many years for the case of finite temperature 
gauge theories, where L = (3 = l/T. The effective potential for the Polyakov loop is easily 
calculated to lowest order in perturbation theory; for a pure gauge theory it is given by 

—2 00 Tr aP™ 



7l 2 (3 



n=l 



where the trace of P in the adjoint representation is given by Tr^P = |Tr^P| 2 — 1. This 
effective potential is minimized when P = zl where z G Z(N), indicating that Z(N) 



symmetry is spontaneously broken at high temperatures where the one-loop expression is 
valid. Note that when P = zl, the effective potential is just the free energy density of a 
blackbody with 2 (iV 2 — 1) degrees of freedom. It can easily be shown that any additional 
fields also act to break Z(N) symmetry at high temperatures, corresponding to small L — (3. 

There are two broad approaches to maintaining Z(N) symmetry for small L. The first 
approach deforms the pure gauge theory by adding additional terms to the gauge action 



[9 
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13|. The general form for such a deformation is 

/oo 
d 3 x Y, a kTr A P{x) k . (13) 
fc=i 

Such terms are often referred to as double-trace deformations. If the coefficients are 
sufficiently large, they will counteract the effects of the one-loop effective potential, and 
Z(N) symmetry will hold for small L. Strictly speaking, only the first [N/2] terms are 
necessary to ensure confinement. It is easy to prove that for a classical Polyakov loop 
P, the conditions TrpP k = with 1 < k < [N/2] determine the unique set of Polyakov 
loop eigenvalues that constitute a confining solution, i.e., one for which Tr^P = for all 
representations with &r ^ Ijj]. The explicit solution is interesting: up to a factor necessary 
to ensure det P=l, the eigenvalues of P are given by the set of iV'th roots of unity, which are 
permuted by a global Z{N) symmetry transformation. The effective potential associated 
with S is given approximately by 



V^ loop (P, 0) = —j £ -±- + £ a k Tr A P k . (14) 

n P n =l n k=l 



For sufficiently large and positive values of the a^'s, the confined phase yields the lowest 
value of Vi-i 00 p. However, a rich phase structure emerges from the minimization of V\-\ oop . 
For iV > 3, the effective potential predicts that one or more phases separate the deconfined 
phase from the confined phase. In the case of SU(3), a single new phase is predicted, and 



has been observed in lattice simulations 



For larger values of N, there is a rich set of 
possible phases, including some where Z(N) breaks down to a proper subgroup Z(p). In 
such phases, particles in the fundamental representation are confined, but bound states of 
N/p such particles are not jl2| . 

Lattice simulations of SU(3) and SU(4) agree with the theoretical predictions based on 
effective potential arguments [9J. The phase diagram of SU(3) as a function of T = L" 1 and 
di has three phases: the confined phase, the deconfined phase, and another phase, the skewed 
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phase, where charge-conjugation symmetry is spontaneously broken. In the case of SU(4), a 
sufficiently large value of a\ leads to a partially-confining phase where Z(4) is spontaneously 
broken to Z(2). Particles with k = 1 are confined in this phase, i.e., (TrpP (x)) = 0, but 
particles with k = 2 are not, as indicated by {TrpP 2 (x)) ^ 0. An important result obtained 
from the lattice simulation of SU(3) is that the small- L confining region, where semiclassical 
methods yield confinement, are smoothly connected to the conventional large- L confining 
region. 

Another approach to preserving Z(N) symmetry for small L uses fermions in the adjoint 



representation with periodic boundary conditions in the compact direction [10] . In this 
case, it is somewhat misleading to use {3 as a synonym for L, because the transfer matrix for 
evolution in the compact direction is not positive-definite. Periodic boundary conditions in 
the timelike direction imply that the generating function of the ensemble, i.e., the partition 
function, is given by 



Tr 



(15) 



where F is the fermion number. This graded ensemble, familiar from supersymmetry, can 
be obtained from an ensemble Tr [exp — /3H)] with chemical potential \x by the re- 
placement — > in. This system can be viewed as a gauge theory with periodic boundary 
conditions in one compact spatial direction of length L — (3, and the transfer matrix in the 
time direction is positive-definite, 

The use of periodic boundary conditions for the adjoint fermions dramatically changes 
their contribution to the Polyakov loop effective potential. In perturbation theory, the 
replacement f3fi — > in shifts the Matsubara frequencies from (3u n = (2n + 1) 7T to (3u n = 2nn. 
The one loop effective potential is now essentially that of a bosonic field, but with an overall 
negative sign due to fermi statistics 15) . The sum of the effective potential for the fermions 
plus that of the gauge bosons gives 



1 ~ Tr A P n 



Vx-ioov (P, m, N f ) = — 4 £ 2N f (3 2 m 2 K 2 (n/3m) - — 



2 



(16) 



^ 4 ^i n 2 

where Nf is the number of adjoint Dirac fermions and m is their mass. Note that the first 
term in brackets, due to the fermions, is positive for every value of n, while the second term, 
due to the gauge bosons, is negative. 

The largest contribution to the effective potential at high temperatures is typically from 
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the n = 1 term, which can be written simply as 

1 



7T 2 /3 4 



2N f f3 2 m 2 K 2 (/3m) - 2 \Tr F P\ 2 - 1 (17) 



where the overall sign depends only on Nf and (3m [16]. If Nf > 1 and (3m is sufficiently 
small, this term will favor TrpP = 0. On the other hand, if (3m is sufficiently large, a value 
of P from the center, Z(N), is preferred. Note that an = 1 super Yang-Mills theory 
would correspond to Nf = 1/2 and m = 0, giving a vanishing perturbative contribution for 
all n 111, [3]. In that case, non-perturbative effects lead to a confining effective potential for 
all values of (3. In the case of Nf > 1, each term in the effective potential will change sign 
in succession as m is lowered towards zero. For larger values of N, this leads to a cascade of 
phases separating the confined and deconfined phases. As m increases, it becomes favorable 
that TrpP n 7^ for successive values of n. If A" is even, the first phase after the confined 
phase will be a phase with Z(N/2) symmetry. As m increases, the last phase before reaching 
the deconfined phase will have Z(2) symmetry, in which k = 1 states are confined, but all 
states with higher k are not. Lattice simulations of 57/(3) with periodic adjoint fermions 
are completely consistent with this picture 19[]. For AT > 3, there are generally phases 
intermediate between the confined and deconfined phases which are not of the partially- 
confined type. Careful numerical analysis appears to be necessary on a case-by-case basis 
to determine the phase structure for each value of A" 20]. 

There are some interesting additional issues arising when periodic adjoint fermions are 
used to obtain Z(N) symmetry for small L. There are strong indications from strong- 



coupling lattice calculations 



21] and the closely related Polyakov-Nambu-Jona Lasinio 



(PNJL) models [22j that the mass m that appears in the effective potential V(P) should 
be regarded as a constituent mass that includes the substantial effects of chiral symmetry 
breaking. In strong-coupling lattice calculations and PNJL models, this effect is responsi- 
ble for the coupling of P and ipip in simulations of QCD at finite temperature. Thus it is 
important that recent simulations of ST/ (3) with Nf = 2 flavors of adjoint fermions show 
explicitly that the Z(A^)-invariant confined phase is regained when the fermion mass is suf- 
ficiently small. In fact, the simulations show the behavior expected on the basis of effective 
potential arguments: a skewed phase separates the confined phase and deconfined phase. 
However, this raises another issue: For a simple double-trace deformation, the small- L and 
large-L regions are smoothly connected. Is this also the case if periodic adjoint fermions 
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are used? A semi-phenomenological analysis based on a PNJL model |23| suggests that the 
answer is yes. However, there may be problems in reaching sufficienly light fermion masses 
in lattice simulations to confirm this. 

In addition to the overall phase structure, semiclassical methods may also be used to 
determine string tensions for small L; the same methods are used with double-trace de- 
formations and with periodic adjoint fermions. There are two classes of string tensions, 
depending on whether the Polyakov loop or the Wilson loop is used. For simplicity, we will 
refer to the former electric and the latter as magnetic, borrowing the terminology 

from finite temperature gauge theory. Electric string tensions are calculable perturbatively 
in the high-temperature confining region from small fluctuations about the confining mini- 



mum of the effective potential 



16 



24| . The m = limit has the simple form 



"PLY 



(2iv f -i) fl 2 r 2 

3iV 



3 csc^ 



7TA 



N 



- 1 



(18) 



and is a good approximation for (3m <C 1. The confining minimum of the effective potential 
breaks SU(N) to [/(l)^ -1 . This remaining unbroken Abelian gauge group naively seems 
to preclude spatial confinement, in the sense of area law behavior for spatial Wilson loops. 
However, as first discussed by Polyakov in the case of an SU (2) Higgs model in (2 + 1)- 
dimensional gauge systems, instantons can lead to nonperturbative confinement j^. There 
are N different fundamental monopoles in the gauge theory on R 3 x S 1 . This behavior was 
first seen in the construction of finite temperature instantons (calorons) with non-trivial 



Polyakov loop behavior at infinity |26H30l] . As in the case of Polyakov's original model, the 
string tension is proportional to exp (— S monopo i e ), where S monopo i e is the action of a single 
monopole. The other factors are essentially determined by the renormalization group up to 



an overall numerical factor 
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3l|. 



There remain, however, a number of interesting issues which require further analytic 
study, including string tension scaling laws for Wilson loops, the mechanism for chiral sym- 
metry breaking, the approach to the large- N limit of SU(N), geometries with more than 
one compact direction 11| . and gauge groups other than SU(N). Lattice simulations have 
the potential to provide much useful information, for example on the connection between 
confinement at small L and large L. Because of the small value of the running coupling 
constant at small L, gauge theories on R 3 x S 1 are a natural area where lattice simulations 
and analytic calculation can work together effectively. 
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IV. CONFINEMENT VIA REAL-SPACE RENORMALIZATION GROUP 



SU(N) gauge theories in four dimensions have a single, Gaussian ultra-violet fixed point 
at g 2 = 0. This means that QCD with massless fermions has no free parameters once a 
renormalization-group invariant length scale A is defined. However, non-perturbative fea- 
tures can so far only be determined in lattice simulations. The natural framework for deal- 
ing with such problems is a Wilsonian RG blocking procedure bridging the different scale 
regimes. A variety of lattice gauge theory approaches have been developed, including tech- 
niques for finding the action along the Wilsonian renormalization group trajectory ("perfect 
action" ) and Monte Carlo renormalization group techniques. An alternative approach which 
has been extensively developed by Tomboulis is based on approximate real-space renormal- 
ization group transformations which satisfy rigorous bounds 32|-l35l|. The general strategy 



is to employ approximate but easily explicitly computable renormalization group transfor- 
mations that provide bounds on the ratios of partition functions, leading to results on free 
energy differences. Sufficiently strong bounds constrain will constrain the exact free energy 
differences. 

The Migdal-Kadanoff approximation provides a starting point for rigorous analysis 
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391 ] . The Migdal-Kadanoff procedure uses a combination of decimation and bond-moving 
to construct an approximate renormalization group which gives a rigorous lower bound on 
free energies. Decimation refers to the procedure of integrating over some of the variable 
of a lattice system, e.g. a lattice spin system such as the Ising model or a lattice quantum 
field theory. Suppose the partition function for such a system is defined on lattice A as 
Z\ = Tr^exp(-S) where S is the Euclidean action for a field theory or f3H for a classical 
system. The Tr\ indicates that all variables on the lattice A are to be integrated over. One 
way to construct a real-space renormalization group transformation is to integrate over some 
of the variables, leaving only the variables on a sub-lattice A' and a new action 5" for those 
variables. This may be written as 

exp(-5") = Tr A/A / exp(-S') (19) 

where the trace is over variables on the complement A/A' of A. A typical choice for A' is to 
go from a lattice with lattice spacing a to a sub-lattice of the same type with lattice spacing 
ba with b > 1. However, this procedure can only be carried analytically for the simplest 
of models, notably d — 1 spin systems and d = 2 gauge theories with local interactions. 
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In higher dimensions, new, complicated quasi-local interactions are generated. The Migdal- 
Kadanoff procedure uses bond-moving to restrict the terms generated by decimation. Bond 
moving refers to modifying the lattice action from S to S + AS in such a way that the 
average action is unchanged: Tr A e s AS = 0. In the simplest case, AS decouples some of 
the variables on A/A' from those on A' to eliminate unwanted couplings. This procedure 
gives a bound Z'a> on the original partition function 

Z' A , = Tr A exp(-S - AS) = Z A {^\ > Z A (20) 

where Jensen's inequality has been used in the form ( e_A5 ) A > e^~ As ^ A = 1. With a 
good choice for AS, this procedure leads to a relatively simple set of renormalization group 
equations for the coupling constants of the model. This procedure can be generalized in 
principle to include blocking of variables, in which a variable on A' represents a weighted 
average of variables on A. The Migdal-Kadanoff approach has been usefull y applie d in 



lattice gauge theory to give an understanding; of renormalization group flows 
the deconfmement transition 



40 



42 



4l] and 



43J in lattice gauge theories. 
Tomboulis has generalized the Migdal-Kadanoff procedure in such a way that not only is 
the partition function left invariant, but other quantities are left invariant as well. Suppose 
the goal is to calculate or estimate an expectation value 

(O) = %^ (21) 

The goal is then to construct a renormalization group procedure consisting of a set of steps 

a ba -> b 2 a ->■ ... -»■ b n a 

A — > A (1) — > — > ... — > 

such that 

(O) = = ... = Za - [ ° ] = ... = Za - [ ° ] (22) 
Z \ Z \ Z k 

After a sufficiently large number of renormalization group transformations, an asymptot- 
ically free theory will reach the strong-coupling region near the infra-red fixed point at 
g 2 = oo, where cluster expansion methods can be used reliably. Note that the construction 
is specific to a given operator O: different operators will give rise to different renormalization 
group flows. In principle, the difference between the two flows would be irrelevant in the 
continuum limit. 
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The prototypical application of this approach is to the twisted partition function of an 
577(2) lattice gauge theory 44] . The twisted partition function Z^ "* on an Li x L 2 x L 3 x L4 



lattice A is obtained from the usual partition function Z\ by modifying the action such that 
plaquettes in the x — y plane with fixed values x$ and y$ of x and y are multiplied in the 
action by —1 for all values of z and t: U xy (x ,y , z,t) —> —U X y(x ,y ,z,t). The twisted 
partition function was originally defined in the continuum by 't Hooft and plays a central 



role in his analysis of the allowed phases of gauge theories |45|, 4(|. The vortex free energy 
Fa is defined by the ratio of the two partition functions: exp [—-Fa] = Z A ~^ /Za; in this ratio 
the bulk free energy term proportional to space-time volume cancels, 't Hooft's analysis 
shows that confinement holds if 

F A ~ cL 3 L 4 e~ pLlL2 (23) 

with p > as the infinite space-time volume limit is taken with L1L2 ^> log (L3L4). For 
lattice gauge theories, a theorem due to Tomboulis and Yaffe proves rigorously that the 
Wilson loop has area law behavior if p is non-zero : 

A/L X L 2 



(W(C)><2[i (i_ e - 



F A 



(24) 



where A is the area associated with C [41 

There remain some unresolved questions with this construction. There has been some 
criticism of both the implicit nature of the renormalization group construction and its pos- 
sible inability to distinguish between asymptotically-free theories like 577(2) and the non- 



asymptotically free U(l) lattice gauge theory 48-[50]. It is true that the renormalization 
group flow is defined implicitly, and there is not yet a complete demonstration or working 
implementation of the specific renormalization group procedure proposed. However, this is 
related to the second, more physical criticism, which is based on the inability of Migdal- 
Kadanoff renormalization group equations to differentiate between non-Abelian and Abelian 
groups in the flow equations. This is in fact a very subtle issue. The four- dimensional £7(1) 
lattice gauge theory has two phases: a weak-coupling phase which is a lattice version of a 
free photon theory, and a strong-coupling phase where magnetic monopoles are responsible 
for confinement. After decades of effort, lattice simulations have not yet conclusively deter- 
mined the order of the transition between these phases. It is believed that the weak-coupling 
phase consists of a line of renormalization group fixed points, similar to that found in the 
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d = 2 XY model, but no reliable renormalization group calculation shows this. The Migdal- 
Kadanoff renormalization group shows a dramatic change in the renormalization group flow 
near the known critical point: the renormalization group beta function becomes very small, 
becoming proportional to exp (— c/g 2 ). This behavior suggests that the renormalization 
group is trying to account for magnetic monopole effects. Exactly the same behavior is seen 
when the Migdal-Kadanoff scheme is applied to the d = 2 XY model J5I]. In this case, there 
is a reliable renormalization group calculation available, but it depends on the introduction 
of a vortex fugacity and a renormalization group flow with the coupling and the vortex fu- 
gacity as independent couplings. In the generalized Migdal-Kadanoff renormalization group 
developed by Tomboulis, the domain of decimation parameters is enlarged relative to the 
original Migdal-Kadaonff scheme. The issue of the existence of a fixed point versus a quasi- 
fixed point which fails by exponentially small corrections to be a true fixed point in the 
U(l) gauge theory thus requires renewed investigation in this more general context before 
it is settled. Thus the original debate points to the large, difficult issue of accounting for 
the infrared effects of topological excitations, particularly within a renormalization group 
framework. 

Beyond issues concerning renormalization group flows in pure gauge theories, there are 
important issues arising from the inclusion of fermions within a real-space renormalization 
group framework. Fermions appear in functional integrals as Grassmann variables, i.e. anti- 
commuting c-numbers, and cannot be simulated or approximated in the same way as bosonic 
variables. Fermion fields that appear only quadratically in the action can be formally in- 
tegrated over, yielding a functional determinant that acts as an additional weight in the 
bosonic functional integral. In lattice simulations, this functional determinant is included, 
essentially as an additional, non-local interaction that is expensive to simulate. In the con- 
text of real-space renormalization group calculations, any partial integration over fermions 
with light masses would introduce non-local interactions difficult to incorporate into a renor- 
malization group calculation. Furthermore, the bounds on partition functions that hold for 
bosonic functional integrals no longer hold for fermionic degrees of freedom. Nevertheless, 
;here has been some progress. Some analytical work has been carried out for free fermions 
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53), and there has been more recent work in connection with the rooting problem of 



staggered lattice fermions [54]. In work reported at the workshop, Tomboulis discussed a 
blocking technique for fermions that is capable of producing a non-trivial fixed point in 
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SU(2) lattice gauge theory for a sufficiently large number of flavors Nf. The existence of 
such a fixed point is indicated by perturbation theory, but the fixed point is generally located 
outside the region of validity of perturbation theory. The presence of a non-trivial infrared 
fixed point indicates the existence of a "conformal window" which is of great interest for 
studies of particle physics beyond the standard model 
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56|. 



V. CONFINEMENT VIA THE FUNCTIONAL RENORMALIZATION GROUP 
AND DYSON-SCHWINGER EQUATIONS 

The closely related functional renormalization group and Dyson- Schwinger approaches 
to quark confinement have much in common with the research on confinement on R 3 x S 1 
described in section 3, including the construction of the Polyakov loop effective potential. 
However, these approaches focus on the infrared properties of conventional gauge and ghost 
propagators. In the functional renormalization grouo approach the renormalization group is 
used directly in attempting to determine the non-perturbative infrared properties of gauge 
theories from the renormalization group flow equation {5?]]. The Schwinger- Dyson approach 
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59[ aims at the calculation of the same quantities by a self-consistent solution of the 
Schwinger- Dyson equations, which are the equations of motion for correlation functions. 
Although beyond the scope of this review, both approaches can be used to study the interplay 



of confinement and chiral symmetry breaking 



59 



60|, with results similar to those obtained 



from more phenomenological approaches such as the PNJL model 



22]. 



The functional renormalization group can be elegantly described for any field theory by 



the RG flow equation for the effective action T [(f)] |6ll . l62j . The effective action T [(f)] is the 
generator of one-particle irreducible (1PI) vertices, where <fi is shorthand notation for one 
or more classical fields, one for each quantum field. The n-point 1PI vertices are obtained 
from T [(f)] by functional differentiation: 

T^(xi,...,x n )= ° (25) 

6(f) (xi) ...5(f) (x n ) 

The 2-point function T^(xi, £2) is the full propagator. The functional renormalization group 
generalizes T [(f)] to [(f)] , which is identical to V [(f)], except that only momenta with q 2 > k 2 
are included in the n-point vertices. This exclusion of infrared modes is accomplished by 
inserting an infrared cutoff function Rk into the functional integral that defines T [(f>], so that 
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by definition T [(f)] = T [(f)]. The RG equation for T k [(f)] is 

£™ = Hl?iSF*M (26) 

where the trace involves an integration over momenta or coordinates as well as summation 
over any internal indices. For the study of spontaneous symmetry breaking, the key quantity 
is typically the effective potential V (0), obtained by setting to a space-time independent 
value, and dividing T [0] by the volume f2 of space-time: V ((f)) = T [0] /Q. 

The functional renormalization group formalism can be applied to gauge theories with 
some modifications of the procedure for scalar fields: a gauge-fixing term must be added to 
the action and ghost fields must be introduced. In an asymptotically free theory, one can 
approximate 1^ by the classical action S for large k, and evolve the functional RG equation 
down to k — 0. The effective potential for A A is given by V (At) = Y [A4] /Q, where A 4 is a 
constant, representing a constant Polyakov loop. This effective action can be evaluated using 
the background field method, and the effective potential is obtained within the functional 
renormalization group approach from 

V k (At) = ^ (27) 

in the limit k — > 0. In principal, the functional renormalization group is exact, but in 
practice, approximations are necessary. A scaling solutions is assumed in the infrared, in 
which the transverse gauge field propagator scales as 1/ (p 2 ) 1+KA and the ghost propagator 
scales as 1/ (p 2 ) 1+Kc . Then at sufficiently low temperature, the effective potential is given 
by 

V (A 4 ) = |^ (1 + k a ) + \ - (1 + «c)| • ^Tr log [-D 2 [A 4 ]} (28) 

where — D 2 [AA is the standard scalar propagator for a particle at finite temperature in a 
background field A4. This should be compared with the standard high-temperature expres- 
sion 

V (A 4 ) = {^-l}>lo g [-D 2 [A,]} . (29) 
If the solution of the functional renormalization group equations satisfies 

^ (1 + ka) + \ ~ (1 + kg) < (30) 

then the sign of the effective potential changes between high- and low-temperature, leading 
to the same confining minimum at low temperature found at high temperature in section 3. 
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This scenario was studied in Braun, Gies, and Pawlowski [57| using Landau-DeWitt gauge, 
and the values ka = —1.19... and = 0.595... were obtained, indicating confinement. 

A closely related approach is based on the self-consistent solution of the Dyson-Schwinger 
equations in the infrared 58] . The Dyson-Schwinger equations are an infinite set of coupled 
equations relating n-point correlation functions to correlation functions of higher order. 
Again using a shorthand notation, the Dyson-Schwinger equations may be derived from 



functional integration by parts 



5 



(x 1 )..4(xn)e- s[lt)] 



(31) 



5<f>(y) 

The simplest case of n = reduces to (6S/8<p(y)) = 0, the statement that the classical 
equation of motion holds on the average. In both the functional renormalization group and 
Dyson-Schwinger approaches, an infinite tower of coupled functional equations for the cor- 
relation functions are solved, and both approaches include the same physics, albeit captured 
in different ways. At present, studies using the Schwinger-Dyson approach have restricted 
themselves mainly to the infrared regime relying on a general power counting analysis of 
the possible IR fixed points of the full system of equations of quenched QCD, i.e., without 
quarks, in the Landau gauge. Applying this approach to gauge theories requires exten- 
sive graphical analysis to identify the most important contributions to the Dyson-Schwinger 
equations in the infrared. Choice of gauge is important in simplifying the analysis. One 
finds both a decoupling solution that does not feature any infrared enhancement and a scal- 
ing solution where infrared divergences are induced in the ghost sector {ssj. In the scaling 
solution, a long range interaction in gauge-dependent local correlation functions in the quark 
sector arises from a strong kinematic infrared divergence of the quark-gluon vertex. In this 
approach, the Wilson loop is then given by an infinite series of these local Green functions 
and the infrared divergence of the leading term, given by the quark-quark scattering kernel, 
yields in the quenched approximation area law behavior [59]. Recently work indicates that 
the corresponding long range interaction is a universal property of the gauge sector 63]. 
Although there are some significant differences at a technical level, the Dyson-Schwinger 
and functional renormalization group approaches are very similar, and in Landau gauge, the 



infrared behavior found for the ghost and gluon propagators is essentially the same 58l. 164 1 . 
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VI. CONCLUSIONS 



The problem of quark confinement is an important and vital subject in theoretical physics, 
connected not just to hadronic physics, but to many areas of physics by concepts and tech- 
niques. It is the quintessential problem of dimensional transmutation, where renormalization 
introduces a length scale into a classicaly scale-invariant system. Fundamental questions of 
mechanism, i.e., the cause of confinement, have not yet been answered. In addition to ques- 
tions of mechanism, there are also questions surrounding the implications of confinements, 
e.g., string tension scaling laws and other string properties. 

Approaches have been developed in the last few years that are shedding new light on these 
problems. The discovery of a rich set of new phases in gauge theories on R 3 x S 1 shows that 
there are new phenomena that remain to be uncovered. The topics discussed here, despite 
differences of approach and emphasis, have more in common with each other than the use of 
the renormalization group. For example, recent work using the functional renormalization 
group and the Dyson-Schwinger equation shares with work on R 3 x S* 1 an emphasis on 
the effective potential for the Polyakov model. As discussed in section 3, there are issues 
concerning magnetic monopoles raised by recent work using real-space renormalization group 
methods that are related to recent work on R 3 x S 1 . Beyond common technical issues, there 
are two basic themes all the approaches discussed here have in common: an incorporation of 
both the physics of continuum gauge theories indicated by the renormalization group and the 
increasingly precise results of lattice gauge theory. This marks a significant departure from 
previous approaches to confinement, which tended to neglect either lattice or continuum 
results, and should be welcomed as genuine progress towards greater understanding. 
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